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Abstract 

The twist-bend nematic phase, IVtb> may be viewed as a heliconical molecular arrangement in 
which the director n precesses uniformly about an extra director field, t. It corresponds to a nematic 
ground state exhibiting nanoscale periodic modulation. To demonstrate the stability of this phase 
from the elastic point of view, a natural extension of the Frank elastic energy density is proposed. 
The elastic energy density is built in terms of the elements of symmetry of the new phase in which 
intervene the components of these director fields together with the usual Cartesian tensors. It 
is shown that the ground state corresponds to a deformed state for which K 22 > A 33 . When the 
elastic free energy is interpreted in analogy with the Landau theory, it is demonstrated the existence 
of a second order phase transition between the usual and the twist-bend nematic phase, driven by a 
new elastic parameter playing a role similar to the one of the main dielectric anisotropy of classical 
nematics and being closely related to the bulk compression modulus representing the pseudo-layers 
of twist-bend nematic phases. A phase transition and the value of the nanoscale pitch are predicted 
in accordance to experimental results. 
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The recently discovered twist-bend nematic phase, JVtb, has been theoretically predicted 


by Meyer [1] and Dozov Q, and has been experimentally evidenced by a number of studies 


15] . By speculating that the bend elastic constant, K33, may become negative, Dozov [2] 


employed a simple fourth-order model for the bulk free energy to predict the existence of 
two different periodic one-dimensional textures in the nematic phase. The first one was 
called splay-bend texture, characterized by a local bend that changes periodically its sign. 
The other one is a continuous conical twist-bend (TB) texture, in which the director n 
rotates along an axis forming a revolution cone with aperture 6. In achiral system, this 
latter texture is twofold degenerated, allowing both right-hand and left-hand twists. By us¬ 
ing techniques of small-angle x-ray scattering, modulated differential scanning calorimetry, 
it was recently concluded that the low-temperature mesophase of CB7CB is a new uniax¬ 
ial nematic phase whose director distribution is composed of twist-bend deformations [3]- 
Transmission electron (TEM) and polarized optical microscopy have been used to experi¬ 
mentally demonstrate the existence of the TB nematic phase 15]. Sctructural observations 
confirming the existence of a TB ground state have been also carried out by TEM, to¬ 
gether with measurements of the director cone angle and the full pitch of the director helix, 
indicating a strong coupling between the molecular and the director bend [J|. These ex¬ 
perimental identifications permitted to describe this nematic order as a true liquid-crystal 
phase with a new type of order. It is then a twist-bend modulated phase formed by achiral 
molecules, in which the director follows an oblique helicoid, maintaining a constant oblique 
angle 0 < 9 < tt/2, with the helix axis 1_5]. If the helix axis is along z, then the director n 
may be written as n = (sin 6 cos (p, sin 6 sin Lp, cos 6), where ip is the azimuthal angle, given 
by ip = qz = (27 r/p)z, where q is the modulus of the wave vector and p is the nanoscale pitch 
of the helix. In this way, this new ground state may be faced as an intermediate structural 
phase between the (usual) uniaxial and the chiral nematic phases 15|. It can be also faced 
as a non-uniform ground state, characterized by the presence of a local spontaneous bend 
and a very low value of the bend elastic constant K33. To account for the elastic properties 
of the phase, Virga 0 proposed an intrinsically quadratic elastic theory to explain Ntb 
phases with an extra director field, in which there are two variants of the helical nematic 
phase, as predicted in Ref. Q, with an helix axis t: one in which the helix winds upward 
and another one in which it winds downward, as shown in Fig. [0 This theory does not 
require a negative bend elastic constant, as suggested by the quartic theory proposed in 
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FIG. 1. One of the variants of the helical nematic phase in which the helix can wind upwards or 
downwards. The vertical axis defines the direction of t. 


Ref. [2] nor the existence of a locally ferroelectric phase, as required by the theory proposed 
in Refs. jl7, 18|, which is also quadratic but deals with an effective bend constant arising 
from the flexoelectric coupling. 

In this letter, a quadratic elastic theory also based on the existence of an extra director, 


;ory alsc 
f. [L6|, is 


represented by the helix axis t, as in Ref. [16], is proposed. Here, however, this extra director 
plays the role of an external held that couples with the nematic director n. Differently from 
what is done in Ref. |l6j, in which the IVtb is assumed as existing, the reference state used 


here is the undeformed one, characterized only by the uniform part of the free energy density. 
The aim is to demonstrate that a phase like the twisted-nematic one may be stabilized, i.e., 
the ground state of the system corresponds to the minimum of a quadratic elastic theory built 
using only the existing elements of symmetry of the phase. To the best of our knowledge, 
this is the first theory capable of predicting the TB phase as a ground state, stable phase. 
It can also predict the TB to regular nematic phase transition and estimate parameters 
connected to experimental results. Nevertheless, both theories reduce to the usual Frank’s 
one for ordinary nematics if the extra director vanishes. 


The starting point is to consider a crystal characterized by a director n and a torsion t 
related to the chiral twisted collective arrangement. One assumes no polar order, in such a 
manner that n is the usual nematic director. The free energy of a crystal of this type, /, 
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depends on the director field n. According to general rules 

/(n) = / 0 (n) + Lijriij + ^K idtk ,in id n Kh 


( 1 ) 


where the dots are terms of higher order in n hJ = drii/dxj, and the summation convention 
has been assumed. The quantity /o(n), independent of can be decomposed, at the 


lowest order, as 


/o( n ) = /i-^(n-t ) 2 + 


( 2 ) 


where f\ is the uniform part of the energy of the usual nematic phase, if one treats the vector 
t as an external field, as proposed above. In this framework, the parameter rj plays the role 


of the main dielectric anisotropy of the nematic director with respect to t 


20]. In this sense, 


?? may be viewed as the strength of the deformation of n about t. It is then connected to 


, 0 ). 


the compression modulus of the pseudo-bulk layers B (as discussed in Ref. 

The tensor of second order has to be decomposed in terms of the elements of symmetry 
of the phase, n and t, of the identity tensor and of the antisymmetric tensor of elements Sij 
and £ijk, respectively. Standard calculations give 

Lij = AiJiiHj + A 2 riitj + Atfinj + A^tjtj + A 5 Sij + A 6 rik£kij- (3) 

Since the medium is globally non-polar, the condition 


/( n ) = /(— n ) ( 4 ) 

holds. It requires that L i3 has to be odd in n, which implies that A\ = A 4 = A$ = 0. 
Furthermore, since |n| = 1, i.e riiTii = 1, it follows that rqrq.., = 0. Consequently, the term 
connected to A 2 does not play any role. The elastic energy density linear in the deformation, 
related to the spontaneous deformation of the phase under consideration, is 

-hp' Tlij tr | tj flj Tlj j T hv2Tlk £kij ^ • (h) 

where the phenomenological constants appearing into Eq. (J3]) have been renamed as = 
—A 3 and k 2 = Aq. A simple calculation shows that 


U rij riij = —t • [n x (V x n)] and rik £kijni,j = n • (V x n), (6) 

and the linear term in the deformation of the elastic energy can be rewritten in covariant 
form as 


Lij n i,j = fiq t • [n x (V x n)] + n 2 n • (V x n). 


( 7 ) 
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One notices that while the term related to n • (V x n) is a chiral term, present in stan¬ 
dard cholesteric liquid crystals, the first term is peculiar to the heliconical phase under 
consideration. In the case in which t = u 2 , the director held n may be expressed as 

n = [cos <p(z) u x + sin <p(z) u y )] sin 9 + cos 9 u,, ( 8 ) 

where Uj,, u y , u 2 are the unit vectors along the Cartesian axes. As stated before, the director 
held given by Eq. (JHD corresponds to a precession of n around z with constant orientation 
of n with respect to t. In this case, one easily shows that 

t ■ [n x (V x n)] = 0 and n • (V x n) = — q sin 2 9 (9) 


are constant quantities. 

The fourth-rank tensor, K t jki, comparing in the elastic energy ([I]), may be decomposed 
by following the standard procedure (see e.g., Ref. [19]), namely: 


A ijkl k^TljfliSik T k^S^jS^l T ^7 T ^S^il^jk T P1 T 

T ^ Q^ititjSkl T tktlfiij ) T Qs^i^k^jl T ~(j4{titl : Sjk T tjtk&il) 

+ q^tjtAk + nUejki. ( 10 ) 

The terms involving Aq yield the usual Frank contribution. The other terms are also decom¬ 
posed in the usual manner and the total elastic energy density may be written as: 

/ = fo ~ ^h( n • t) 2 + «i t • [n x (V x n)] + k 2 n • (V x n) + ^AT U (V • n) 2 + ^K 22 [n • (V x n)] 2 

+ i/l 33 (n x V x n) 2 - (K 22 + K 2 4 )V ■ (nV • n + n x V x n) 

+ Pi[t • (n x V x n)] 2 + gi[t ■ V(t • n)] 2 + q 2 [ t • V(n • t)(V • n)] 

+ g 3 [V(t • n)] 2 + <74[(t ■ V)n] 2 + g 5 [V(n • t) ■ (t • V)n] - riV(n • t) • (V x n). 


For the physical state represented by the director given by Eq. fl 8 lh the terms connected 
with gi, q 2 , g 3 , < 75 , and rq do not contribute, because n • t = cos#, which is independent of 
z. The only term that survives is the one connected with qq, because: 

[(t • V)n] 2 = q 2 sin 2 9. 

For what concerns the usual Frank terms, the non-vanishing ones are 


( 11 ) 


n x V x n) = q 2 sin 2 9 cos 2 9 and (n ■ V x n) 2 = q 2 sin 4 6 
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After collecting all the non-vanishing terms, one obtains for the elastic energy density the 
following expression: 


f(q, x) = fi- ^?7(1 -x)- n 2 qx + i I< 22 q 2 x 2 

+ ^K 33 q 2 x(l -x) + q 4 q 2 x, ( 12 ) 

where x = sin 2 9. The procedure is now to minimize f(q,x ) with respect to x and q. By 
imposing df /dx = 0 and df /dq = 0 , one obtains 


A 33 + 2q 4 — k 2 \J (A 33 + 2 ^ 4 Jjrj 

A 22 — K'i'i 


for the cone angle of n with t, and 


(13) 


= Vray < 14 > 

for the modulus of the wave vector. The director profile associate to x* and q* corresponds 
to a minimum of the free energy density of an unlimited nematic sample only if 

'a 2 /' 


dx 2 


> 0 , 


x 


and the Hessian determinant defined by 


, * (d 2 fd 2 f d 2 f d 2 f \ 

^ ,q ) “ {jh^-dq 2 ~ > ' 


Simple calculations yield 


dx 2 


A 22 — A; i3 

= q— --— > 0 


x* ,q* 


and 


A 33 + 2 g 4 

H(x*, q*) = , y i/iA".,; + 2g 4 ) - )/(A'., + 2g 4 ) > 0. 


(15) 


(16) 


(17) 


(18) 


These are the requirements to be fulfilled in order to have an equilibrium state. From 
Eqs. (ED and m, one concludes that two scenarios are possible. Since q* is real, one 
obtains that a stable ground state may exist if rj < 0 and K 33 + 2 q A < 0. In this case, 
J](K 22 — K 33 ) / (K 33 + 2g 4 ) > 0 when K 22 > K 33 . On the other hand, another stable ground 
state may be obtained if r] > 0 and K 33 + 2q 4: > 0. Again, rj(K 22 — K 33 )/(K 33 + 2 q 4 ) > 0 when 
K 22 > K 33 . Thus, both situations are physically possible, implying that the twist elastic 
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FIG. 2. Elastic energy density F = f — f\ (in arbitrary units), given by Eq. IfTQlb as a function of 
x in three situations: 1) for A > 0 and B > 0 (dotted line), 2) A = 0 and B > 0 (dashed line), and 
3) A < 0 and B > 0 (solid line). For illustrative purposes, the curves have been drawn for q > 0. 


constant is such that K 2 2 > A' 33 , as observed experimentally (see, e.g.,Ref. 22]). Now, the 
analysis can be cast in the framework of a simple Landau-like expansion of the free energy 
if one rewrites Eq. (fl2l) in a still more compact form as 


f{x) = fi ~ | + Ax + Bx 2 , (19) 

in which A = —F 2 + 77/2 + K 33 /2 and B = A/2, with the following compact notation 
introduced: F 2 = n 2 q, K 33 = (A " 33 + 2 q 4 )q 2 , and A = ( K 22 — K 33 )q 2 . Simple calculations 
show that / is minimum for x = —A/(2B), as given by (fT3l) when (fT4l) is taken into account. 
This corresponds to a stable state if B > 0, i.e., A > 0 or Ii 22 > K 33 . The twist bend 
phase, which corresponds to x 7 ^ 0, is energetically favored only if A < 0 (see Fig. [2]). Thus, 
there is a phase transition from the usual nematic phase (x = 0 ) to a twist-bend nematic 
x 7 ^ 0 for A = 0, i.e., when 2 n 2 q — rj + (K 33 + 2q 4 )q*, which, by using (fT4lh becomes 
rj — r) c — k 2 2 /{K 33 + 2q 4 ). For this critical value, the wavevector of the modulated phase 
assumes the value q c = |K 2 /(A 33 + 2g 4 )|. In addition, A = rj — s/rfiy. may be used to drive the 
transition according to the values of rj. Since the elastic energy density is being interpreted 
in the framework of a Landau theory, one can conclude that a second order phase transition 
occurs when ij = r] c , if x = sin 2 6 is used as the “order parameter” for the IVtb phase. This 
critical value of r] separates the usual from the twist-bend nematic phase. In this way as 
well, the parameter r] may play a role similar to the temperature such that A > 0, for 
ij > rj c , and A < 0, for q < q c . The order of magnitude of q in the pure IVtb phase may 
be easily estimated if, for simplicity, one neglects the cholesteric contribution (k 2 = 0), 
since the phase is not chiral. From Eq. (fT3|) . one obtains q = q 2 K 22 sin 2 9 — K 33 q 2 sin 2 9 or 
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7] = B / sin 2 9 — ( q 2 K 33/sin 2 9) sin 4 9, where B = K 22 q 2 sin 4 6 is the compression modulus 
mentioned above 0 - Since K 22 > /I33, one can assume that 7 7 ~ B/sm 2 9. By using 
B « 2 kPa, if, from (TT41) . we consider that q = yjB/{K 22 — AT 33 ) sin 2 9 , with 9 = 16 °, K 22 ~ 
10 -12 N, and Jl 33 w 10 ~ 13 N, one obtains that p 0 = 2n/q « 10 rim, in good agreement with 


experimental results for helix pitch of the twist-bend phase 4 , [ 15 J, 2lJ. These estimations 


indicate that the new parameter not only can be expressed as a combination of known 
parameters characterizing the IVtb phase, but is also found to have an order of magnitude 
in good agreement with the available experimental results. 


Some concluding remarks are in order. The experimentally discovered twist-bend nematic 
phase may be described from the elastic point of view by means of a quadratic theory, which 
can be constructed as the usual one, i.e., by considering only the elements of symmetry 
enough to characterize the system. For the twist-bend phase, these elements are the director, 
n, and the helix vector t, the latter playing a geometrical role similar to the one of an external 
field. In addition to the usual Frank constants, the elastic theory requires the existence of 
three parameters, namely, rj, which plays a role similar to the one of the main dielectric 
anisotropy in usual nematics, k 2 , which is the coefficient of the twist distortion, peculiar 
to the heliconical phase, and q 4 , a parameter that renormalizes the elastic constant A'33. 
The model can be interpreted in analogy with the Landau phenomenological approach for 
phase transitions. This analogy permits to treat the cone angle 9 as the main ingredient 
for an "order parameter" defined by x = sin 2 9 . The twist-nematic phase is then the one 
for which x 7^ 0, and is energetically favored with respect to the usual nematic phase 
(x = 0 ) when q = rj c . There is a critical value of the parameter r) pointing a second order 
transition between the two stable phases allowed by symmetry considerations. This scenario 
is characterized by the exigence that K 22 > /I33 always, as experimentally evidenced. In 
addition, the possibility of having negative values of /i 33 , as argued by some authors, may 
not be excluded but it is not mandatory. Indeed, one of the possible stable phases may 
exist for rj < 0 and /I33 + 2 q± < 0, but this condition may be satisfied for /I33 > 0 when 
the elastic parameter q 4 is negative enough, i.e., when the term q 4 [( t • V)n] 2 = q 4 q 2 sin 2 d is 
energetically favored, which is a physically sound possibility. 
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